Extracting the ground-state spin of a quantum dot from the conductance peaks in a 

parallel magnetic field at finite temperature 
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We derive a closed expression for the finite-temperature conductance of a Coulomb-blockade 
quantum dot in the presence of an exchange interaction and a parallel magnetic field. Parallel- 
field dependence of Coulomb-blockade peak position has been used to determine experimentally the 
ground-state spin of quantum dots. We find that for a realistic value of the exchange interaction, 
the peak motion can be significantly affected at temperatures as low as kT ~ 0.1 A, with A being 
the mean level spacing in the dot. This temperature effect can lead to misidentification of the 
ground-state spin when a level crossing occurs at low fields. We propose an improved method to 
determine unambiguously the ground-state spin. This method takes into account level crossings and 
temperature effects at a finite exchange interaction. 

PACS numbers: 73.23.Hk, 73.63 Kv, 71.70.Ej, 73.40 Gk 



The ground-state (GS) spin of a mesoscopic structure 
such as a quantum dot has been a topic of major inter- 
est Q. In the absence of interactions, the total kinetic 
energy plus confining potential energy of the electrons is 
minimized in eigenstates with minimal spin. However, 
the exchange interaction is minimized in eigenstates of 
maximal spin. The GS spin is therefore determined by 
the balance between these two competing energy terms. 

The GS spin of a mesoscopic structure is difficult to 
measure. In devices in which charging energy is impor- 
tant, Coulomb blockade (CB) serves as a useful tool to ex- 
plore GS spin P|. In particular, the dependence of a CB 
peak position on a parallel magnetic field has been used 
to determine the GS spin in quantum dots j^, , 
metallic nanoparticles 7] and carbon nanotubes jal- At 
temperatures kT that are much smaller than the single- 
particle mean level spacing A of the dot, the conductance 
peak position is determined by the change in the GS en- 
ergy as an electron is added to the dot. The Zeeman 
coupling to an in-plane magnetic field B\\ is described 
by gi^BB\\S\\, where g is the gyromagnetic factor (taken 
to be positive), /xs is the Bohr magneton, and is the 
component of the total spin along the field. Thus, in 
the presence of a parallel field, the GS energy of the dot 
with GS spin S is shifted down by an amount qubB^S, 
corresponding to the state of lowest magnetic quantum 
number M = —S. The CB peak position acquires a shift 
— (7/zb-B|| A5, where A^* is the change in the GS spin as 
the number of electrons increases from to -I- 1 . The 
CB peak positions are then expected to exhibit a cer- 
tain pattern of ^gus/'^ slopes versus describing the 
increase or decrease of the GS spin by 1/2 0. 

An excited state in a dot with energy E* (at zero mag- 
netic field) and spin S* > S will acquire a larger Zeeman 
energy shift and will cross the GS level Eqs at a value of 
the parallel field given by giXBB\\{S* - S) ^ E* - Eqs- 
If the magnetic field is further increased, the excited en- 



ergy level will become the new GS for the dot, resulting 
in an abrupt change in the CB peak position slope. The 
new slope is still determined by the change of the dot's 
total spin but it now involves the spin of the correspond- 
ing zero-field excited state instead of the GS spin. In the 
presence of an exchange interaction, higher spin states 
shift down in energy, and crossings are expected to occur 
at smaller values of -B|| . 

Slopes of =Fg/is/2 and traces of crossings are clearly 
seen in the small Si dots of Ref. for small values of 



of Ref. 



A < 0.3). However, for the weakly coupled dots 
the observed slopes of the peak spacings at 
small are essentially flat (see Fig. 9), and were difficult 
to interpret. Similar flat slopes were seen in Ref. Q. 

Here we investigate CB peak motion in quantum dots 
under an applied parallel magnetic field and at finite tem- 
perature. Remarkably, we find that signatures of level 
crossings can be almost completely washed out at tem- 
peratures as low as kT ^ 0.1 A. In general, we observe 
that the peak positions become flat at sufficiently low 
(for which the Zeeman energy g^BB\\ is below kT) |^. 
However, when a crossing occurs at small values of B\\, 
the peak position remains flat up to the crossing point. 
These findings suggest a possible explanation of the re- 
suhs of Refs. 0] and 0, for which kT ~ 0.15 A and 
kT ^ 0.2 A, respectively. If the slopes beyond the flat 
section are used, the GS total spin can be misidentified. 
We propose an improved method in which the CB peak 
position and peak height data are both fitted to a two- 
transition model that includes an excited state in the 
dot with either A^ or A^ -I- 1 electrons. We show that this 
method can detect all relevant crossings at temperatures 
kT ~ 0.1 — 0.15A, allowing for a correct assignment of 
GS spins in a systematic way. 

We assume a quantum dot that is weakly coupled to 
leads, such that F ^ kT and F ^ A, where F is a typical 
tunneling width of an electron from the leads into the 
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dot. In this limit we can use a rate-equations approach 
to calculate the linear conductance G in the pres ence of 
residual interactions ^3 • The results of Ref. |lOl | can be 
generalized to include a parallel magnetic field. 

A closed solution for the conductance is not always pos- 
sible. An important case for which an explicit solution 
exists is the universal Hamiltonian, obtained for a chaotic 
or diffusive dot in the limit of a large Thouless conduc- 
tance 0, ^1 . In the presence of an applied magnetic 
field parallel to the plane of the dot, the Hamilto- 
nian is ^ = X^A ^-^"-^ + e^n^/2C - JsS^ -I- g^Bi?||S||, 
where e\ are spin-degenerate single-particle levels, C 
is the dot's capacitance, Jg is the exchange interac- 
tion constant, and S is the total spin of the dot. The 
term g^BB\\S\\ describes the Zeeman coupling to the 
field -B|| , where 5*11 is the total spin projection along the 
field direction. The orbital-level occupations spin 
S and spin projection M along the field direction are 
all good quantum numbers, and the eigenenergies are 



^%]m = Ex CA^A + e^Ny2C - JsS{S +1)+ gtisB^^M. 

Here we generalize the expression for the conductance 
G in the vicinity of an ^ + 1 CB peak, obtained in 
Ref. to the case of an applied parallel field. We find 

A 

where T is the average width of a level, 
.gA= 2f-ir\r5^/(ri^ -hT'^) are the single-level (di- 

mensionless) conductances, and are the single-level 
tunneling widths to decay from level A to the left 
(right) lead. The thermal weights w'f^ and w'"^^ collect 
the contributions to the conductance from processes 
in which an electron tunnels into an empty or singly 
occupied level A, respectively. Unlike the case without 
the Zeeman term, explicit summations over the magnetic 
quantum numbers remain, and we find 
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Here Pn,s.m = e"'^[^"-s+^".s.Afl/Z is the probabil- 
ity of finding the dot with A^ electrons, spin S and 
spin projection M, with Fjq^s being the free energy 
of A^ non-interacting electrons with total spin 5*, and 
Un^sm = e'N^/2G - J,5(5 + 1) + gfiBB^M - ipN. 
The Fermi-Dirac distribution f{(,ss'cr) evaluated at 
Cs,S',. = ^A- Js {S'{S' + 1)- S{S + l))-eF-|-gMsS||^/2 
and a/2 = M' — M is the spin projection of the elec- 
tron that tunnels into the dot. ep = ep + e.GgVg/G 
is an effective Fermi energy where Vg is the gate volt- 
age and Cg the dot-gate capacitance. The coeffi- 
cients 6a,jv,s = {^{fix- l){nx - 2))^^^ and cx^n+i,S' = 
(^"'a("'A ~ 5, are spin- and particle-projected 

quantities whose explicit expressions are given in Ref. 

Using a realization of the single-particle Hamiltonian, 
we calculate the conductance from Eqs. Q, ^ and as 
a function of ep, and find its maximum to determine the 
peak position. In Fig. ^a) we show calculated peak po- 
sition versus g^BB\\/ IS. for four consecutive peaks (A, B, 
C, D) and for a realistic exchange constant of Js = 0.3 A 
At kT — 0.01 A (dashed fines) the expected slopes 
of igiJ-B /2 are observed. Assuming a starting value of 
S = below peak A and using the observed peak position 
slopes, we assign GS spin values of 5 = 0, 1/2, 1, 1/2, as 
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FIG. 1: (a) CB peak position versus gfisB^^/A for ,7$ = 0.3 A 
and four consecutive peaks at kT = 0.01 A (dashed lines), 
and kT = 0.1 A (solid lines), (b) CB peak spacing versus 
(;/i_B_B||/A for the peaks shown in (a). The dotted lines are 
the expected slopes of igfiB and 0. 

electrons are added to the dot. At the higher tempera- 
ture of kT = 0.1 A (solid lines) the peak position curves 
become flat near = (i.e., for g^BB\i < kT) but the 
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FIG. 2: (a) CB peak position versus (;/i_BS||/A for Jg = 0.3 A 
and for 5 consecutive peaks at temperatures of kT — 0.01 A 
(dashed lines), kT = 0.1 A (dash-dotted lines), and kT = 
0.15 A (solid lines), (b) CB peak spacing versus qubB^^/A for 
the peaks shown in (a). The dotted lines describe slopes of 
±g^B- Inset: the probability distributions for — > 1, 1/2 
3/2 and 1^2 crossings to occur at parallel field By and for 
Js = 0.3 A. 



low-temperature slopes can still be identified. 

Fig. nib) shows the corresponding peak spacings ver- 
sus 5/iBi?||/A. Peak positions with alternating slopes 
lead to peak spacings with slopes of ±g^B- However, two 
parallel consecutive peak positions result in a zero-slope 
peak spacing (e.g., B — A) and indicate that the GS spin 
increases (or decreases) twice in a row, e.g., — > 1/2 ^ 1. 
We see that the peak-spacing slopes do not change much 
at the higher temperature of kT = O.IA and thus GS 
spins can be correctly inferred at this temperature. 

Fig. [Sfa) shows peak position versus gfisB^^/A for a 
different set of five consecutive peaks (E, F, G, H, I) for 
Js = 0.3 A and temperatures of kT = 0.01 A (dashed 
lines), kT = 0.1 A (dash-dotted lines) and kT = 0.15 A 
(solid lines). At kT = 0.01 A, level crossings are clearly 
observed in the form of kinks in the peak position versus 
For example, the dot between peaks H and I has 
5 = in its zero-field GS, but at g^BB\\/A ^ 0.3, an 
5 = 1 excited state crosses the S — state and becomes 
the new GS. The slope of the peak position versus is 
still given by gfiB^S, but now AS* describes the change 
of the spin involving the new GS of the dot in which 
the level crossing occurs. Thus the 5/1^/2 slope of peak 



position H (1/2 — » transition) changes into a —gfis/'^ 
slope (1/2 — ^ 1 transition) following the level crossing at 
g^BB\\/A ~ 0.3. Signatures of level crossings are also 
seen in the peak spacings of Fig.Efb) at kT = 0.01 A. 

However, at the higher (but still low) temperatures of 
kT = 0.1 A and kT = 0.15 A, signatures of level crossings 
have almost completely disappeared. We observe that 
when a level crossing occurs at a low value of B\\ , the peak 
position and peak spacing at kT = 0.15 A are essentially 
flat up to the crossing point. Flat slopes were observed 
at low fields for weakly coupled dots in Refs. 11 and [| 
(see Fig. 9 and Fig. 5, respectively) and are likely the 
result of temperature effects and/or level crossings. If 
the peak position and peak spacing are flat up to the 
crossing point, it is not possible to identify the slopes 
that are necessary for the correct GS spin identification. 
At higher values of the field, the peak spacing acquires a 
slope but this slope reflects the GS spin of the dot after 
the crossing. If such slopes are used to determine the 
zero-field GS spin, the spin value can be misassigned. For 
example, the dot between peaks H and I can be assigned 
a spin of 5 = 1 instead of its true GS spin of 5 = 0. 

Spin misassignmcnt is likely when a level crossing oc- 
curs at relatively low values of the field, reflecting a zero- 
field low-lying excited state with spin S* > S. The prob- 
ability distributions for various level crossings to occur 
at a field B\^ are shown in the inset of Fig. El^b) for 
Js — 0.3 A. In particular, the probability of a — + 1 
spin crossing at low _B|| is quite substantial (^ 0.3 for 
(7^B-B||/A < 0.3). Thus the occurrences of 5 = 1 ground 
states are likely to be overestimated. 

To avoid spin misassignmcnt, we propose an improved 
method to determine the GS spin of the dot from par- 
allel field measurements at temperatures kT < 0.15 A. 
The method is based on the observation that at these 
temperatures the conductance is dominated by the con- 
tribution from two groups of transitions between the iV- 
and {N + l)-electron dots, obtained by considering one 
level with spin S in one of the dots and two levels with 
spins 5*0 = S — 1/2 and S*! = 5 + 1/2 in the other dot 
(e.g., 1/2 —> and 1/2 —> 1 transitions). In this case, 
the rate equations can be solved explicitly for any resid- 
ual interactions Assuming the relevant transitions 
from the N- to the {N + l)-electron dots are S ^ Sq 
and S ^ Si with reduced tunneling widths Fq and Fi, 
respectively, we find 
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spin S in the A''-electron dot, and is the probabiHty 
to find the dot in its TV-electron state with spin S and 
spin projection M. A similar expression can be derived 
when the relevant N ^ N + 1 transitions are So ^ S 
and 5*1 S. This two-transition model takes into ac- 
count the possibility of level crossing (versus -B|| ) in the 
N- or N + 1-electron dot. The conductance formula 
is valid for arbitrary residual interactions, but is a low- 
temperature approximation kT < 0.15 A. On the other 
hand, the conductance given by Eqs. (^), |2Jl and ||2Jl is 
only valid for the universal Hamiltonian but is exact for 
arbitrarily large temperatures. 
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FIG. 3: Results of the two-transition model Eq. (|1J (symbols) 
are fitted to the full calculation data of Fig. |21 (solid lines) for 
both peak positions (panel (a)) and peak heights (panel (b)). 
Simultaneous fitting of both peak positions and peak heights 
is needed to properly extract the values of the total spin from 
the data. Circles indicate crossing positions. 

To determine unambiguously the GS spin, we fit simul- 
taneously the measured peak position and peak height 
data versus parallel field using the above two-transition 
model. In a given range of _B|| , we assume a certain 
spin S (for either the N~ or the (TV -I- l)-electron dot) 
and fit four parameters: the two conductances Gi cx 
r'r^/(r' -l-F^) and the two many-body energy differences 



(e.g.. 



.(JV+l) (N) 



Sg if S* is the spin of the A^-electron dot). 



We validate our method using as data sets the kT — 
0.15 A peak positions of Fig. I^^a) and their correspond- 
ing peak heights, shown by the solid lines in Fig. OJa) 
andlSfb), respectively. To emulate the experimental sit- 
uation, we treat the spins values and the crossing posi- 
tions as unknown. The fit of these data sets to the two- 
transition model (symbols in Fig.O is found to be good. 
In particular, the corresponding spin values and crossing 
positions (dashed circles in Fig. (SJ^a) and (b)) can be ex- 
tracted systematically at a temperature of kT = 0.15 A. 



We have also verified that the parameters we extract 
from the fit agree with their input values in the universal 
Hamiltonian. While it is possible to fit separately the 
peak position and peak height data, the best-fit param- 
eters are often found to be unphysical. 

Our method can also be used to extract the lowest 
energy level (at zero field) for a given spin value that is 
higher than the GS spin, as long as this level becomes 
the GS of the dot at some finite parallel field. 

In conclusion, we have derived a closed expression for 
the conductance in a parallel magnetic field of a CB quan- 
tum dot that is described by the universal Hamiltonian. 
For typical low temperatures of kT > 0.1 A used in the 
GaAs quantum dots, signatures of level crossings in a par- 
allel field can be almost completely washed out. When 
such a crossing occurs at low values of the field, the GS 
spin of the dot, as determined from the peak position 
motion in a parallel field, can be misidentified. We have 
described an improved method in which the measured 
peak position and peak height data are fitted to a finite- 
temperature two-transition model that takes into account 
level crossings. The model is valid for any interaction in 
the dot, and thus could be generally used as to extract 
the GS spin in experimental situations where the lowest 
attainable temperature is ~ 0.1 A. 
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